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The propagation of magnetoelastic waves in a magnetic insulator having a 
nonuniform internal magnetic field is examined in the geometrical optics 
approximation. Hamilton's ray path equations are obtained from the 
slowness relation for the medium, and it is shown that for YIG there is a 
substantial focusing action in the rod configuration commonly used for 
magnetic delay line experiments. When external field shaping is used to 
produce a minimum internal field at the midpoint of the rod it is found 
that divergence of the magnetoelastic waves is to be expected. 

I. INTRODUCTION 

In a number of experiments, 1 - 23 propagation of magnetoelastic waves 
has been observed in discs and rods of yttrium iron garnet. Coupling is 
provided through an internal field variation along the direction of 
propagation, radially in a disc and axially in a rod. This permits excita- 
tion of the wave in a region of small wave vector, 4,6 where the magnetic 
field can couple to the magnetization, with subsequent tapering into the 
magnetoelastic crossover region. The demagnetizing field also varies in 
magnitude and in orientation across the direction of propagation. In 
regions where the wave vector is large it is appropriate to consider the 
effects of this field inhomogeneity in terms of geometrical optics, and it 
is to be expected that refraction of the magnetoelastic waves will occur. 

II. THE SLOWNESS RELATION AND GROUP VELOCITY 

In a cubic crystal with a dc magnetic field applied along a [100] axis 
x% and, for simplicity, assumed elastically isotropic propagation of mag- 
netoelastic waves is governed by the set of equations 6 
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iuM Xl + u Hk M X3 — iybk cos 8 R t > = 
(w fffc + o) M sin 2 6)M Xl - iwM X2 - iybk (cos 20 R t + sin 20 Bi) = 

(co 2 - c,V)# t - i(bk/ P M) cos 20 Jlf Xl = (1) 
(co 2 - *¥)£i< - i(bk/ P M) cos 6M X2 = 
(to 2 - C|V)fii - i(bk/ P M) sin 20 M Xl = 0, 
where 

co fl , = 7 (H + # e *aV) 
COAf = t47tM. 

The wave vector k is assumed to lie in a (100) plane at an angle 6 
with the dc field. M X1 , M X2 are transverse components of the magnetic 
moment referred to axes along [100] directions, and B t > , Rt , Ri are 
transverse and longitudinal components of elastic displacement. The 
saturation magnetization is denoted by M and the mass density by p. 
Transverse and longitudinal elastic wave velocities are represented by 
c t and ci respectively, and b is the second magnetoelastic constant, 
generally designated byb 2 .H is the internal dc magnetic field, H ex is the 
exchange field, and a the lattice constant. In what follows it will be 
assumed that the crystal has sufficient magnetoelastic isotropy (6i « b 2 ) 
that (1 ) is valid for a magnetic field applied at a small angle to the [100] 
axis and for propagation in any azimuthal direction. 

Upon elimination of variables in (1), the secular equation is found to 
be 

o*' 1 a\ ( 2 2 ^ ( &fc ) 

i2(co,/c,0) = (to, — co J 



" MHt , , (2) 

f to 2 cos 2 6 , i»H k cos 2 26 . uh** sin 2 2d \ _ Q 

'{utf - u* + to, 2 - co 2 "*" to, 2 - to 2 j 



where 



j, 2 = wj/ t (w» t + com sin 0) 



cow = coj — 



, (6A;) 2 2 



P MH, 

2 . 2j.2 

COj 

CO J 



COS 



2 = c t 2 k* 



2 = e,V. 
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This equation relates the wave vector k or the wave slowness vector 
k/co = k/kv p h to co and 0, a relation which may be displayed graphically 
by a dispersion diagram (see Fig. 1 ). In coordinates fa ,k 2 ,fa the slow- 
ness relation (2) defines a "wave vector" 7 or "slowness" 8 surface for 
each value of co. Since the magnetoelastic dispersion curves (see Fig. 1 ) 
have four branches and the dispersion relation is independent of azi- 
muthal angle, the wave vector surface is a surface of revolution about a- 3 
and comprises four sheets. For example, a vertical section through the 
sheet of the wave vector surface corresponding to branch III appears, 
at co tt co w , as shown in Fig. 2. The group velocity vector 



V„ = V fc co = - 



dQ/dco 



(3) 



is proportional to the gradient of 12 and is therefore normal to the wave 
vector surface, as shown in Fig. 2, the sense of the normal being de- 
termined by the requirement that the angle between V and k be less 
than x/2. This means that except in the special cases = or t/2, the 
group velocity vector is not exactly parallel to the wave vector; and a 
wave packet does not move in a direction normal to its phase fronts, a 
phenomenon which is characteristic of anisotropic media. 



w = C t K 




Fig. 1 — Magnetoelastic dispersion diagram. 
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PURE ELASTIC WAVES 

MAGNETOELASTIC WAVES 

Fig. 2 — Wave vector surface for branch III of the dispersion diagram, at 



III. THE EIKONAL EQUATION AND THE RAY EQUATIONS 

It is assumed that the magnetic field varies in both magnitude and 
direction from point to point in the medium, but is sufficiently strong at 
all points to saturate the magnetization. The geometrical optics ap- 
proximation is appropriate when the magnetic field is almost constant 
over regions comparable with a wavelength in dimension, so that a solu- 
tion to the magnetoelastic equations having the form 

M(r)«* w 

R(r)e* w 

appears over a small region as a plane wave with relatively slowly vary- 
ing amplitude. If the assumed solutions are introduced into the equations 
of motion and spatial derivatives of R and M are neglected, equations 
(1) are obtained with | V^ | substituted for \k \, and the angle between 
the "local" wave vector V^ and the local magnetic field is substituted 
for 0. With the same substitutions, the slowness relation (2 ) reduces to a 
first-order partial differential equation for the phase function \f/, 
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Sl(u,p> , Xt) = 

*-*. -1,2,3. (4) 

0.li 

This is the eikonal equation or equation of geometrical optics. In the 
neighborhood of a singularity of the medium the approximations made 
in deriving the eikonal equation sometimes break down. For the case 
of electromagnetic propagation in a ferrite, Seidel has shown that 
singularities of this kind occur because of the appearance of logarithmic 
derivate coefficients in the field equation. It is not clear whether similar 
singularities exist for the magnetoelastic equations, and no attempt will 
be made to justify rigorously the use here of the geometrical optics ap- 
proximation. 

The standard method of solving (4) is by means of the characteristic 
or ray equations 11 

dXi/dw = dSl/dpi (5a) 

dpi/dw = -dtt/dXi, (5b) 

where w is a parameter. For any set of initial values of p,- , .x\ satisfying 
(4) these equations, which form the basis of Hamiltonian optics, ' de- 
fine a unique curve in the space .r, , x 2 , .r 3 . The significance of this 
curve becomes clear when the equivalence of p, to the component A", 
of the "local" wave vector is recalled. This shows [from (3)] that the 
tangent, 

{dxi:dx 2 :dx 3 ) = h-:7-:7- 

\dpi dpi dp 3 

to any curve defined by (4) and a set of initial conditions is always col- 
linear with the group velocity vector. Therefore the curve, or ray path, 
obtained by integrating (5) describes the trajectory of a wave packet 
launched at a specified point .r, with a specified "local" wave vector 
hi = Pi . The value of the phase function \f/ at any point on the ray path 
is obtained implicitly from 

/ = Pi t- + P2 — + P3 r— , (6) 

dw dpi dp 2 dp 3 

where [from (5a)] dw is related to the increment in ray path length ds 
through the relation 

-<(£) ,+ e) ,+ (fc)7 fc 
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When the wave vector surface has more than one sheet there are 
several initial group velocities corresponding to the same initial wave 
vector direction. These are distinguished by the magnitude of the initial 
wave vector, and a wave packet will therefore trace out different ray 
paths according to the magnitude of the initial wave vector, each path 
corresponding in a local sense to propagation in a mode associated with a 
particular branch of the dispersion diagram. When there is only a slow 
spatial variation of the magnetic field there will be little coupling between 
modes and the different ray paths will maintain their distinct identities. 

The present discussion will be concerned with ray paths corresponding 
to branches I and III of the dispersion diagram. In this case Schlomann 
has shown that, in the region where the uncoupled magnetic and trans- 
verse elastic waves cross, the slowness relation can be written ap- 
proximately as 

(<o - «.)(« - Ctk) - (<T/2)co cr f(e) = 0, 

where a = yb 2 /aM, co cr is the crossover frequency for the uncoupled 
waves, a is the elastic stiffness c« , and 

f(0) = i[(2 - 5 sin 2 + 4 sin 4 6) (1 4- \ umW* sin 4 0) 1 

+ £ WAfwr 1 sin 4 0(3-4 sin 2 0)]- 

At the lower microwave frequencies it can be shown that the k depend- 
ence of oj« due to exchange has a very much smaller effect on the slope 
of the dispersion curves than does the magnetoelastic coupling. If 
exchange is neglected 



= U s = COff 



( 1 + — sin 1 , 



in which w H — yH and the eikonal equation takes the form 

S!( W> P , , .U) = ( Pl 2 + p? + P 3<) ! + £¥^-f=°~ (7) 

ZCt 03 — Cdg Ct 

where is the angle between the vector (pi:p 2 'P3) and the local magnetic 
field. 

In the following section attention will be directed toward rotationally 
symmetric systems, with rays travelling in meridian planes. It is ap- 
propriate, then, to use a cylindrical coordinate system, and (5) and (7), 
which are written in Cartesian coordinates, must be transformed. Since 
the <p component of V^ is zero, (7) becomes 
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^,p r ,p ! ,r,z) = (p r +p 2 ) +--_-- - = (g) 

= v ~ % 

where rj = tan -1 p r /p z and £ is the polar angle of the local dc magnetic 
field (see Fig. 3). In a rotationally symmetric system the ray equations 
(5) transform into 

dr/dw = dSl/dpr, dz/dw = dSl/dp t (9a) 

dp r /dw = -dtl/dr, dpjdw = -dSl/dz (9b) 

IV. PARAXIAL RAY EQUATIONS AND REFRACTION IN CONVERGING AND 
DIVERGING MAGNETIC FIELDS 

The discussion will now be restricted to the paraxial case; that is, only 
rays lying close to the symmetry axis and traveling almost parallel to it 
will be considered. Then 

Furthermore, the rotationally symmetric magnetic field will be assumed 
to be almost parallel to the axis (£ « 1). Since « 1, 

f(6) « 1 - 2.5 2 ; 

and, when coa//co w < 10, 

co a ~ oiii + cj. tf (6 /2). 

If co £3 co g the denominator of the second term in (8) is small and 

8 a co,/(fl) ^ a cojXfl) doi, dd 
dpi 2c i u) — co„ 2c t (co — co,) 2 dd dpi 



d_ 
d.v, 



d_ a_ Us f(e) o_ u.f{8) /cho, dH da. dd\ 

?.c, 2c, co — co, ~ 2c t (co — co„) 2 \dH dXi dd dXi/ ' 

where p t = p r , p z and .r, = r,z. Then 

(* ~ *) 

™=Pl+^±Pi 1 (10a) 

dp r Pz 2C t pz (CO — W B )" 

55.-1, dob) 

d Pi 

where only terms linear in p r /pz and £ have been retained, in accordance 
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Fig. 3 — Orientation of the "local" wave vector vY< relative to the local mag- 
netic held. 

with the paraxial approximation, and w g za wh has been replaced by a> 
in the numerator of the second term in (10a). Similarly 



d£l (TO 



dr 



dr 2c t (u — u) H ) 2 
From (9), (10) and (11) the paraxial ray equations are 



dr 1 / . trojuu ( 



Vr/V* ~ S) \ 

(« — u B y ) 



and 



dH 

tip r _ _a& dr 

dz 2c t (oj — o) fl ) 2 



(ID 



(12) 



(13) 



(14) 



In the paraxial approximation p, is obtained directly from (8), 

Vz ~ c t V ~ 2(« - u> H )) ' 

Consider now the case of a composite magnetic rod, the middle and 
outer sections having saturation magnetizations M and M respectively, 
which is magnetized along its axis (see Fig. 4). The potential function 
for the dipolar field on the center line of the middle section, assuming 

H n 



M' 



2d 



fT 



z M 



M' 



Z=-l Z = l 

Fig. 4 — Composite rod configuration. 
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uniform magnetization, is 

V D (z) - 2tt(M - M')[\ (I - zf + a 2 }* - { (I + zf + a 2 }' + 2z] 

if the end effects of the outer sections are assumed to be negligible. Close 
to the axis the potential function is 14 

M D (r,z) = V D (z) ~ T ~ V D (z). 

Assuming (a/ 1) 2 « 1, this leads to an internal field in the central region 
of the middle section 

HttH - MM - M') %(*- 

, = H r „ MM - M') a_ zr 

* H z ~ i/n Z 2 / 2 



Z 2 V/ 2 2Z 2 , 

(15) 



where only terms up to second order in z/l, r/l have been retained. Equa- 
tion (15) shows that the internal field diverges with increasing | z \ 
when M < M and converges when M > M . This result has been de- 
rived under the assumption of uniform magnetization. Actually the 
magnetization in the rod will itself be nonuniform, and nonuniformity of 
the field will be greater than is shown in (15). 

A plane magnetoelastic wave is assumed to be propagating in the -\-z 
direction at the midpoint of the rod, 2 = 0. Since p r is then zero at this 
point and £ = from (15), it follows from (12) that dr/dz = 0. This 
means that the ray paths are parallel to the axis at z = 0. Elimination 
of p r and p s from (12), (13) and (14) leads to a second-order differential 
equation with variable coefficients for r(z), and the ray path trajectories 
are obtained by solving this equation, subject to the assumed initial 
conditions. In this case a numerical integration is required for a com- 
plete description of the ray paths. If only the direction of refraction is 
required, the following simpler procedure may be used. 

Substitution of (15) into (13) leads to 

r 
dp r 3<tco / N a I 2 

dz Ac t I- 



O 2 2\ 2 

W — C0ff m + - (o>Af — COM') - -J 



where 

com - cow = «y4x(ilf - M') 

(„ . Z(M - M')aV 

"«,„ = T 1 «n + 1 y 
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Over a small range of z close to z = the value of r will not change ap- 
preciably along a ray path and dp T /dz may be integrated directly, giving 

3troj (w M — u) M ') a zr , _. 

fr '"*(.-*)'» (I6) 

to second order in z/l. This shows that the "local" wave vector deflects 
toward the axis with increasing z when M < M and away from the axis 
when M' > M. The corresponding slope of the ray path is found by 
substituting p z and p r from (14) and (16) into (12). That is 

dr _ . j ™ M co H } a zr 

—7- — — A \J t / \~r 



dz 



2(« - «,.)' (l - ^7 \) 

\ 2(oj — 0)ff m )/ 



W//„ 



-- (17) 



up to terms of second order in z/l and r/l, where 



A = 



\ 2(.w — u Hm )/ 



Equation (17) shows that when the field diverges (M < M) mag- 
netoelastic ray paths which are axial at z = will converge, and vice 
versa. This is easily understood in terms of simple physical concepts. 
When the internal field decreases with increasing | z \ it increases with 
increasing r, as shown by (15). For a fixed frequency and propagation 
angle it is clear from Fig. 1 that the "refractive index" kc t /co decreases 
as H increases. If the anisotropy of the dispersion relation is ignored for 
the moment, this means that off-axis rays curve toward the region of 
higher "refractive index" closer to the axis. This isotropic effect is 
enhanced by anisotropy in the dispersion relation. When the wave vector 
is deflected from the magnetic field direction the ray path (defined by 
the group velocity) is deflected even further, as shown in Fig. 2, leading 
to an increased bending of the ray path. 

This enhancement of the refractive effect by anisotropy is represented 
in (17) by the second and third terms under the bracket. In order to 
estimate the magnitude of these effects consider a YIG rod with a/l = 
0.1 and u> Ho = 2ir X 10 9 . For YIG a = 4 X 10 6 sec -1 and a> M = 3.08 X 
10 10 . According to Schlomann 6 one-half the minimum frequency separa- 
tion of the transverse magnetoclastic branches is 

Wmin = (™ cr /2) J « 0W2) 1 = 1-12 X 10 8 . 

If this value is assumed for w — co//,„ the approximations used in obtain- 
ing (16) and (17) are valid when z < 0.1Z, r < 0.01/. At z = 0.1Z, r = 
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50.5 



0.01/ the anisotropic terms in (17) are found to be an order of magnitude 
larger than the isotropic term, and the slope of the ray path is 

t - -7.7 X 10-. 

dz 

On the basis of an extrapolation at this slope, the ray path should inter- 
sect the axis at z ~ 101. The actual intersection would be closer than 
this because the ray path slope changes continuously with z. When the 
signal frequency is shifted closer to w» m , an increased refraction results. 
For example, if 



u> — «» 



= Wmin/3 



the phase velocity of the magnetoelastic wave is, from (8), still within 
a few per cent of the acoustic velocity; but the slope of the ray path at 
z = 0.1Z, r = 0.01Z is now 



dr 

dz 



= -3.2 X 10" 



and the extrapolated intersection point occurs at z ~ 0.3/. This large 
change in refractive power with decreasing o> — co Hm is due to the reso- 
nance denominators in (17) and is an indication of the steep slopes of the 
wave vector surface, Fig. 2, in the vicinity of = 0. The approximations 
used in obtaining (17) are, of course, not valid at resonance but are still 
at least marginally valid in the case considered here. 

V. CONCLUSIONS 

It has been shown that in a uniformly magnetized medium the phase 
and group velocities of a magnetoelastic wave are not collinear except 
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Fig. 5 — Beam steering by means of an auxiliary field. 
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Fig. fi — Examples of lens configurations. 

when the wave vector is either parallel or normal to the magnetic field. 
This effect might be used for steering or switching an ultrasonic beam by 
means of an auxiliary field (see Fig. 5). Since the direction of the wave 
vector remains constant, the phase fronts remain parallel to the trans- 
ducer faces. 

Substantial refraction effects have been shown theoretically to occur 
in a nonuniformly magnetized medium. For the case of a magnetized 
rod it is found that paraxial magnetoelastic rays at frequency co ~ u„ 
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Fig. 7 — Wave vector surface for branch III of the dispersion diagram, at 
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converge when the dipolar field and the applied field are opposing and 
diverge when the fields are aiding. In arriving at these results, the ef- 
fects of losses and scattering due to imperfections have been ignored. By 
a similar analysis it can be shown that an annular permanent magnet or 
a circular coil encircling the rod will act as a converging lens if its field 
aids the applied field and as a diverging lens if the fields are opposing 
(Fig. 6). Paraplanar ray equations can be derived for radial propagation 
in an axially magnetized thin disk at a frequency u ^ !co H (w w + Wn)}*. 
In this case the anisotropic refraction effect is found to oppose the iso- 
tropic effect and can even cause the net refraction to change sign. The 
physical reason for this can be seen by examining the wave vector sur- 
face for this case; see Fig. 7. This shows that a deflection of the wave 
vector away from 6 = r/2 produces a deflection of the group velocity in 
the opposite direction. 
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